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Abstract. This review is focused on the application of specific fluctuation rela- 
tions, such as the Gallavotti-Cohen relation, to ratchet models of a molecular 
motor. A special emphasis is placed on two-states models such as the flashing 
ratchet model. We derive the Gallavotti-Cohen fluctuation relation for these 
models and we discuss some of its implications. 



1. Introduction 

The macroscopic observables of a system at mechanical and thermal equilibrium 
do not vary with time and can be characterized by a finite number of state vari- 
ables. Thermodynamics imposes a priori constraints on the average values of these 
state variables that are satisfied regardless of the specific nature of the system. 
Because of this property of time-invariance, equilibrium is often imagined as be- 
ing associated with stillness and frozen dynamics. This, of course, is not true: 
a system, even at thermodynamic equilibrium, is constantly evolving from one 
micro-configuration to another. This endless motion at the microscopic level can 
be probed macroscopically by measuring the fluctuations of some physical observ- 
ables, the most famous example being Brownian Motion. Equilibrium fluctuations 
are perfectly well explained by the classical laws of statistical mechanics. 

Brownian motion, its nature, its origins, have been a puzzle to physicists 
during the XlXth century pQ. One paradoxical issue was the following: can one 
rectify this random fluctuating motion and use it to perform some work ? If such a 
rectifying device could be constructed then work would be extracted from a single 
heat reservoir, contradicting the second law of thermodynamics. The most famous 
example of a mechanical system that may play the role of such a Maxwell's demon 
is the ratchet and pawl system, presented by Feynman in Chapter 46 of the first 
volume of his Physics Lectures [2]. A related paradox was proposed in 1950 by 
Brillouin J3J: consider an electrical circuit composed of a diode and a resistor at 
temperature T . The current in the circuit has a zero average value, but because of 
thermal noise, it exhibits non-vanishing fluctuations (known as Johnson-Nyquist 
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noise). Can the diode be used to rectify the current, allowing us to use it to perform 
work? The solutions of these paradoxes is now well-known: thermal fluctuations 
are a universal phenomenon and all systems at a given temperature are subject to 
it. The rectifying device, whatever it is, is also subject to Brownian motion and 
undergoes some unavoidable fluctuations. If the signal to be rectified is produced 
at temperature T and the rectifier is at the same temperature then the thermal 
fluctuations of the rectifier render it totally ineffective and the second law is saved. 
If the rectifier is at a lower temperature then this apparatus can indeed generate 
work: but now there are two heat sources at different temperatures, in accordance 
with the second law. 

In recent years, a renewed interest has arisen in ratchet models in the context 
of non-equilibrium statistical physics. Many ratchet models exist (for a review see 
[3]). In one kind of ratchet called brownian motors, an association of two parti- 
cles, one asymmetric and another one symmetric, can rectify thermal fluctuations 
provided that the two particles are in contact with heat baths at different tempera- 
tures [S]. In another family of ratchet models, the need of two heat baths is removed 
by coupling the ratchet to some external 'agent' (e.g., a chemical reaction) that 
continuously drives the system out of equilibrium. In this case too, under certain 
conditions, work can be extracted (6[[7l[8]. Again, there is no contradiction with 
thermodynamics here: the system is far from equilibrium and the ratchet plays 
simply the role of a transducer between the energy put in by the agent (e.g. chem- 
ical energy) and the mechanical work extracted. The analysis of the energetics of 
such devices far from equilibrium requires concepts that go beyond the classical 
laws of thermodynamics and this remains a very challenging and important open 
issue JSQTJ. 

Biophysics provides numerous examples of systems far from equilibrium. For 
example, a significant part of the eukaryotic cellular traffic relies on 'motor' pro- 
teins that move along filaments similar in function to railway tracks or freeways 
(kinesins and dyneins move along tubulin filaments; myosins move along actin fil- 
aments) [TT|. The filaments are periodic (of period ~ 10nm) and have a fairly rigid 
structure; they are also polar: a given motor always moves in the same direction. 
These molecular motors appear in a variety of biological contexts: muscular con- 
traction, cell division, cellular traffic, material transport along the axons of nerve 
cells... A biological cell forms a crowded environment in which molecular motors 
work together and with other proteins. In these conditions, collective effects arise 
due to interactions between motors [12]. In many cases, these interactions can 
be modeled as excluded volume interactions, and for this reason, the behavior of 
an ensemble of motors in a low dimension space can be described by dynamical 
models similar to the ones developed for traffic problems [T31 U3] ■ In the following, 
we focus on single molecular motor properties, in order to clarify in this simpler 
situation, the dynamics and the energetics of this system far from equilibrium. 

Recently, a general organizing principle for non-equilibrium systems has 
emerged which is known under the name of fluctuation relations [THl [TBI 113 ESI HH1 
[20] HJJ. These relations, hold for non-equilibrium steady states but arbitrarily far 
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from equilibrium [22j [23] [24j [25] , they can be seen as macroscopic consequences 
of the invariance under time reversal of the dynamics at the microscopic scale 
[26] . It is interesting to apply these concepts to small systems which can either be 
mechanically driven as biopolymers [27] or chemically driven as enzymes |28) . 

Molecular motors are enzymes which operate stochastically at the level of a 
few molecules, and for this reason they typically undergo large thermal fluctua- 
tions. Generically, single molecular motors have been described theoretically either 
by continuous ratchet models (see e.g. review by Jiilicher et al. [29 ) or by models 
based on master equations on a discrete space [30l EI]. It is possible to give a 
thermodynamic interpretation of the elementary chemical reactions, which occur 
in a discrete and stochastic way in a molecular motor [32] . Such a thermodynamic 
interpretation of chemical transitions has similarities with the thermodynamic 
interpretation of the Langevin equation at the single trajectory level |33) . At a 
macroscopic scale, constraints arise on the operation of these molecular motors, as 
a result of single reaction events occurring stochastically at the microscopic scale. 
These constraints take the form of a fluctuation relation [MJ|32 ( 35, 361 |3"T1 155] , 

The aim of this review is to explain how recent theoretical results in non- 
equilibrium statistical mechanics, namely these fluctuation relations, provide a 
way to understand the non-equilibrium energetics of molecular motors. Note that 
the same framework apply to both biological molecular motors or artificially made 
nanomachines. 

In the first part of this review, we present the two theoretical approaches 
of molecular motors mentioned above, namely the flashing ratchet model and the 
approach based on master equations on a discrete space. In the second part of the 
review, we derive the fluctuation relations for these specific models and discuss 
some consequences of these relations for molecular motors. 



2. Stochastic models of molecular motors 

Molecular motors are enzymes capable of converting chemical energy derived from 
the hydrolysis of adenosine triphosphate (ATP) into mechanical work. There is a 
large diversity of molecular motors, and correspondingly a large number of pro- 
cesses accomplished by these motors within a cell. There are linear motors such as 
kinesins, dyneins, myosins or the RNA polymerases, and rotating motors such as 
the Fi-ATPase motor or the bacterial flagellar motor. These motors drive not only 
intracellular movements, they are also key players in the motility of the cell itself. 
Although, traditionally, these machines were subjects of investigation in biology 
and biochemistry, increasing use of the concepts and techniques of physics in recent 
years have contributed to a quantitative understanding of the fundamental princi- 
ples of operation of these motors. The possibility of exploiting these principles for 
the design of artificial nanomachines has opened up a new field in nanotechnology. 
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2.1. The basic principle 

On the theoretical side, molecular motors have been described by ratchet devices, 
which are systems able to extract useful work out of unbiased random fluctuations 
[55] . A generic model of such a ratchet device is shown in Figure [T] where the 
motor is represented by a small particle that moves in a one-dimensional space. 
At the initial time t — 0, the motor is trapped in one of the wells of a periodic 
asymmetric potential of period a. Between time and tf, the asymmetric potential 
is erased and the particle diffuses freely and isotropically at temperature T. At 
the switching time tf, the asymmetric potential is re- impressed, the motor slides 
down to the nearest potential valley and, because of damping, is trapped in one 
of the wells. The motor has maximal chance to end up in the same well where 
it was at time t = 0. However, it has a small probability to be trapped in the 
well located to the right and, because of the asymmetry of the potential, an even 
smaller probability to end up in the left well. Indeed, in order to be trapped in the 
right well after time tf, the particle must have diffused between t — and t = tf 
over a distance larger than A towards the right. However, to end up in the left well 
it has to diffuse (towards the left) a distance larger than B, which is much less 
probable because B > A. In other words, because the potential is asymmetric, the 
motor has more chances to slide down towards the right: this leads on average to 
a net total current. The particle has used thermal noise to overcome the potential 
barrier and thanks to the asymmetry of the potential it has moved in a well-defined 
direction. 

In order to move the motor consumes r ATP fuel molecules per unit time, 
which are hydrolyzed to ADP + P (see Fig. [5]): 



It is the chemical energy released by ATP-hydrolysis that allows the motor to 
detach itself from the filament it was bound to. This detachment process corre- 
sponds in the basic mechanism to erasing the potential whereas re-attachment of 
the motor at a the switching time tf corresponds to re-impressing the potential. 
Hence the motor undergoes chemistry-driven changes between strongly and weakly 
bound states (attachments and detachments). It is the coupling between chemistry 
and the interaction with the filament that allows the motor to overcome energy 
barriers. Besides, because of the polarity of the filament, the interaction potential 
is asymmetric, allowing directed motion to set in. 

In general, the motor is subject to an external force / ex t which tilts the 
potential. Besides, when ATP is in excess, the chemical potential difference of 
the reaction of ATP hydrolysis, AJ1 = /xatp — Madp — MP becomes positive. More 
precisely, we denote AJ1 = fc^TA/i, where A/i is the normalized chemical potential 
and 



where [..] denotes concentration under experimental conditions and [..] eg denotes 
equilibrium concentrations. 



ATP ^ ADP + P . 




(1) 
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Figure 1. The principle of a Brownian ratchet: by inscribing and 
erasing periodically an asymmetric potential, a directed motion 
of the particle is induced. In this example, the potential is a saw- 
tooth function of period a = A + B. Since B > A, the potential is 
asymmetric. The relative sizes of the probabilities of ending in one 
of the wells are represented by the sizes of the disks in the lowest 
picture. The right and left probabilities being different, this leads 
on average to a net total current J. 



ATP = ADP + P 




Figure 2. Schematic representation of a molecular motor: by 
hydrolyzing ATP, the motor proceeds along the polar filament 
and carries a 'cargo' molecule, which typically exerts a force / on 
the motor. 



A basic problem is then to determine the velocity of the motor v(/ ox t, A/i) 
(mechanical current) and the ATP consumption rate r(/ ex t, A/x) (chemical current) 
as functions of the external mechanical and chemical loads. 
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To summarize, molecular motors move by using the ratchet effect, providing 
an example of a rectification process of Brownian motion. The two basic require- 
ments for obtaining directed motion are: 

(i) an external energy source, provided by the chemical reaction of ATP 
hydrolysis. During this process ATP is consumed and ADP is produced. This 
reaction therefore breaks time-reversal invariance (more technically, it breaks the 
detailed balance condition which holds at equilibrium, and introduces a bias in 
the dynamics of the motor). 

(ii) The polarity of the filament which breaks spatial left-right symmetry and 
allows the motor to move in a well-defined direction. 



2.2. The flashing ratchet model 

In the simplified discussion above, we did not specify the characteristics of the 
switching time tf. Different models are possible: one can consider a deterministi- 
cally forced ratchet in which the binding potential of the motor switches period- 
ically (and even smoothly) from strong-binding to weak-binding. Another possi- 
bility, that we shall discuss here in detail is the flashing ratchet model [29j [4] in 
which the switching of potentials is sudden and occurs at random times generated 
by a Poisson Process. Then, we shall show how to construct a discrete version of 
flashing ratchet model [30l EI] . 

In the flashing ratchet model [221 HH1 SI] , the state of the motor is described 
by a continuous position variable x and by discrete internal states i = 1, 2 corre- 
sponding to different chemical states of the motor. For instance, one could associate 
one state with a configuration where one motor head is bound to the filament (the 
high energy state) and the other state to a configuration where both heads are 
bound (the low energy state). The motor evolves in two time- independent peri- 
odic potentials Ui(x), with i = 1,2. Note that in the basic mechanism discussed 
in the previous section, U\{x) is a saw-tooth potential and U<2(x) is taken to be 
zero. But one can consider the general case where both U\ and U2 are non-zero 
asymmetric potentials of arbitrary shape with a common period a. In figure [31 we 
represent the often studied situation in which XJ\ and U2 are identical saw-tooth 
potentials but slightly shifted with respect to each other along the x axis. 

The dynamics of the motor can be represented by a Langevin equation 

x = -7F - 7 J2 u ^ x )km + VA^W ( 2 ) 

i=l,2 

where is a normalized white-noise and a dichotomous noise that can exist 
in two states 1 and 2. The switching-rates of are position dependent and are 
given by wi(x) (transition from 1 to 2) and u>2(x) (transition from 2 to 1). The 
friction coefficient 7 satisfies the Einstein relation Dq = kBT/^ and F represents 
an external force acting on the motor. The function S^n),i is a Kronecker delta. 
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The probability density for the motor to be at position x at time t and in 
state i is denoted by Pi(x,t), which obeys the following equations 

an n j 

^Qf + ^T = -ui(x)P 1 +o J2 (x)P 2 (3) 

-Qf + -Q^ = M^)Pi-^(x)P 2 , (4) 

where 0J\(x) and u>2(z) are space dependent transition rates, and the local currents 
Ji are defined by: 

with Dq the diffusion coefficient of the motor and F a non-conservative force acting 
on the motor. 

The transition rates can be modeled using standard kinetics for chemical 
reactions applied to each chemical pathway between the two states of the motor 

m- 

wi(ar) = [uj{x) + i,{x)e^]e (lh{x) ' lx)/kBT 1 

uj 2 (x) = [Lo{x)+^{x)}e [U2{x) ~ fx)/kBT , (6) 

where / = Fa/ksT is the normalized force acting on the motor. It is assumed 
that the rates can be decomposed into a contribution proportional to u>(x), which 
is associated with thermal transitions, and a contribution proportional to ip(x) 
corresponding to transitions induced by ATP hydrolysis. Note that the functions 
U)(x) and ip(x) have to be periodic functions but they are otherwise unspecified. 
The form of the rates in the absence of hydrolysis (i.e. when ijj(x) = 0) is chosen 
to satisfy the detailed balance condition 

uji(x) \ k B T J 

The form of the rates associated with the transitions induced by ATP hydrolysis 
(i.e. when oj(x) = 0) is chosen to satisfy a generalized detailed balance condition, 
which is generalized in the sense that it accounts for the exchange of chemical 
energy [351 HE] • In this case, this leads to the condition 

uJi(x) \ kB± J 

Note that the way the force enters the rates is unambiguous in continuous models 
as compared to discrete models, in which the force dependant rates must contain 
unknown load distribution factors [371 133 • One could easily extend the model to 
introduce more chemical pathways |40] or more internal states; such extensions 
are possible but they have not been considered here since they are not essential 
for the present argument. 

When F = and A/i = 0, the system is in equilibrium since the detailed 
balance condition (J7J holds. In this case, the steady state probabilities obey 
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Boltzmann distribution, the currents Ji vanish and there is no average displace- 
ment of the motor. When F and A/x are not simultaneously zero, the detailed 
balance condition ([7J is broken, the system is out of equilibrium and currents are 
present. 

2.3. A discrete ratchet model 

From the continuous model, a simplified effective discrete model can be constructed 
as shown schematically in figure [3] following the procedure outlined in )42j . We as- 




® ® '© — ©- 



Figure 3. The top two curves represent the two time indepen- 
dent periodic potentials U\(x) and Uz{x) of the flashing ratchet 
model. At any position x, vertical transitions are possible between 
the two internal states with rates ijJ\(x) and LU2(x). Below is rep- 
resented the corresponding discrete model, which is obtained by 
considering effective transitions between the minimum of Ui(x) 
(state a) to the minimum of the other potential (state b), with 
rates as shown in the lower part of the figure. 

sume that the motor has a vanishingly small residence time in all states which are 
not minima of the potentials Ui(x) or U^ix), i.e. the time taken to slide down to- 
wards a well is negligible. Switching between the potentials, represented as dashed 
lines in the figure, occurs at finite rates wi(ir) and u>2(x), but only between states 
that are at the same location x. Since downward sliding occurs instantly, the ob- 
servable transitions are effectively from one minimum of one potential (state a) to 
the other minimum of the other potential (state b) . In this way, a discrete hopping 
model on a ID lattice is constructed in which transitions are allowed between even 
and odd sites called a and b. The dynamics of the motor on the linear discrete 
lattice is as follows: the motor hops from one site to neighboring sites, either con- 
suming or producing ATP (see Fig. [4]). The position of the motor is denoted by 
x = nd, where 2d « 8nm is the step size of a kinesin. The even sites (denoted 
by a) are the low-energy state of the motor, whereas the odd sites (denoted by b) 
are its high-energy state; their energy difference is A.E = ksTe, where kg is the 
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Figure 4. A schematic of the evolution of the motor in a plane 
(n, y), where n represents the position of the motor on the filament 
and y is the number of ATP molecules consumed. The even and 
odd sublattices are denoted by a and b, respectively. Note that 
the lattices of a and b sites extend infinitely in both directions 
along the n and y axis. The possible transitions are represented 
with arrows on a particular section of the lattice. 



Boltzmann constant and T is the temperature. Because of the periodicity of the 
filament, all the even (a) sites and all the odd (b) sites are equivalent. The dy- 
namics of the motor is governed by a master equation for the probability, P n (y, t), 
that the motor has consumed y units of ATP and is at site n at time t: 

d t P n {y,t) = -(t7 n + ^„)P„(y,t) 

+ E [^ l n+ xPn+i{y-ht) + ^ l n -xPn-i{y -i,t)], (9) 

i=-l,0,l 

where tj„ = J^i an d = J2i ~^L- Denoted by to^ and ot^ are the transition 
rates for the motor to jump from site n to n — 1 or to n + 1, respectively, with 
/ (= —1,0,1) ATP molecules consumed. 

As we show below, this discrete stochastic model contains the essential fea- 
tures of the original ratchet model while being more amenable to precise mathe- 
matical analysis 43, 44][45l|3T]. In this sense, the discrete model may be regarded 
as a minimal ratchet model. 

2.4. Application of the model to experiments 

2.4.1. Modes of operation and efficiency: From the master equation Eq. [9l one 
can obtain the average velocity of the motor v and its average ATP consumption 
rate, f. One finds explicitly that 

n ~U? a ~ufb — tu a tj(, 

V = 2=r =r ^= ^ (10) 
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Figure 5. Kinesin velocity vs. ATP concentration under an ex- 
ternal force [36]. The solid curves are the fits of our model 
to data from Ref. [46 . From top to down, the plots are for 
F e = —1.05, —3.59, and — 5.63pN, respectively. Inset: Kinesin ve- 
locity vs. force under a fixed ATP concentration. The solid curves 
are fits to the data of Ref. [?B] . From top to down, the plots are 
for [ATP] = 2 mM and 5 fiM. 
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Figure 6. The four modes of operation of a molecular motor 
(such as kinesin) are delimited by v = and r = 0. The lines are 
generated with parameters that we have extracted by fitting the 
data for kinesin in Ref. [46j to our model, and this fit is shown in 

fig. 151 

r (g + gj) tf b + tj b )- (fcj 1 (ga + kg) (n) 

By modelling the dependence of the rates on the force and on the chemical 
potential in a way similar to what was done in Eq. ^ for the flashing ratchet 
model, one obtains a theoretical prediction for the dependence of the velocity and 
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average ATP consumption rate on the force or the ATP concentration, that can be 
compared to experiments. Despite its simplicity, this discrete model can account 
quantitatively for such measurements as shown in figure (J5]) in the case of a kinesin. 

From the two currents v and r, a diagram of operation of the motor (see figure 
[5]) can be constructed, which summarizes the possible thermodynamic modes of 
operation of the motor . This diagram is similar to that given in Ref . , except 
that the present diagram extends to the regime far from equilibrium rather than 
being limited to the linear response regime. Whenever, fv<0 work is performed 
by the motor, whenever rA/i < chemical energy is generated. The motor can 
work in eight different regimes. Four of them are passive and correspond to the 
white regions in Fig. [51 in which there is no energy output from the system, since 
fv>0 and rA/i > 0. □ The four remaining regimes are more interesting since fv 
and rA/i are not of the same sign, which means that some form of transduction 
occurs between the mechanical and chemical forms of energy. More precisely: 

• In Region A of the diagram, where rA/i > and fv < 0, the motor uses the 
chemical energy of ATP to perform mechanical work. This can be understood 
by considering a point on the y-axis of Fig. [6] with A/i > 0. There we expect 
that the motor drifts to the right with v > 0. Now in the presence of a small 
load / < on the motor, we expect that the motor is still going in the same 
direction although the drift is uphill and thus work is performed by the motor 
at a rate W = — fv > 0. This holds as long as / is smaller than the stalling 
force, which defines the other boundary of region A. 

• Similarly, in Region B, where rA/i < and fv > 0, the motor produces ATP 
already in excess from mechanical work. 

• In Region C, where rA/i > and fv < 0, the motor uses ADP in excess to 
perform mechanical work. 

• In Region D, where rA/i < and fv > 0, the motor produces ADP already 
in excess from mechanical work. 

It is interesting to note that the large asymmetry between regions A and C in 
Fig. [6] reflects the fact that kinesin is a unidirectional motor. The diagram also 
illustrates the fact that under usual conditions with A/i ~ 10 — 25, a kinesin uses 
the chemical energy of ATP hydrolysis to produce mechanical work (region A of 
the figure), rather than operating in the other way to synthesize ATP (region B 
of the figure) . 

It is also possible to analyze the mechanical efficiency of this motor, defined 
as the ratio between the mechanical work delivered by the motor divided by the 
chemical energy supplied [37) . Other definitions of efficiency have been consid- 
ered in the literature (such as the efficiency at maximum power [47] or the Stokes 
efficiency) but the advantage of this definition is that it holds arbitrary far from 



1 The case where fv<0 and rA/i < is forbidden by the second law of thermodynamics, because 
it would lead to a negative entropy production. For this reason, there is no point in Fig. l6l which 
corresponds to this case. 
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Figure 7. Curves of equal efficiency r\ within region A (which is 
delimited by the solid line and by the y axis) . The parameters are 
those used in fig. [6] and obtained from the fit of fig. OH From the 
outside to the inside the curves correspond to r/ = 0.2, 77 = 0.3, 
77 = 0.4, rj = 0.5 and 77 = 0.58. The absolute maximum efficiency 
for these parameters is about 59% and is located at A/j ~ 14 and 
/ ~ -4.9. 

equilibrium and it corresponds near to equilibrium to the definition used tradition- 
ally with heat engines. The kinesin operates most efficiently in a range of values 
of A/5 which corresponds well to the typical free energy delivered by the reaction 
of ATP hydrolysis (physiological conditions correspond to A/x ~ 10 — 25). The 
maximum of efficiency is obtained around a single isolated point in the coordi- 
nates (/, A/x) (rather than on a line as in the near equilibrium regime for instance) 
as shown in figure [7J This suggests that kinesin is in fact optimized to operate 
under a load corresponding to a normalized force of about —4.9 in the conditions 
of Fig. [JJ The maximum of efficiency is around 40 — 60%, much higher than the 
typical efficiency in the near equilibrium regime (of the order of 0.03%). The value 
of the maximum efficiency of 40 — 60% agrees well with recent measurements for 
kinesin. 

2.4.2. Violation of Onsager and Einstein relations: Away from equilibrium, we 
expect that Onsager and Einstein relations are no longer valid. To quantify their 
violations, we have introduced in Ref. [36] AA = A12 — A21 to quantify the viola- 
tion of Onsager relations and four "temperature" -like quantities, = Dij/Xij to 
quantify the violation of Einstein relations. All these quantities are defined using 
linear response theory in the vicinity of a non-equilibrium steady state rather than 
near an equilibrium state. Of course, these effective temperatures are not thermo- 
dynamic temperatures: they are merely one of many possible ways to quantify 
deviations of Einstein relations. These and AA are shown in figure |8] as func- 
tions of A/x for the particular case of / C 1 within region A. We observe that 
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FIGURE 8. Effective temperatures Tn (dot-dashed), T 2 \ (dotted), 
T 2 2 (dashed), and AA (solid) vs. A/i in Region A of Fig. [5] for a 
small /. Note that T,-j characterizes the fluctuation- response ratios 
(see text), while AA quantifies the breaking of Onsager symmetry. 



all the Tij start off at Tij — 1 near equilibrium where AA = as expected from 
Onsager relations, whereas for large A[i, T22 diverges exponentially while T21, Tn 
and AA approach constant values [36] . 

2.4.3. Beyond two states model of molecular motors. The discrete two states 
model presented above describes many features of experiments on a single ki- 
nesin such as the average velocity versus force or versus ATP concentration, or 
the average ATP consumption rate. It does not describe however equally well the 
fluctuations of these quantities. As shown in [3D], at least four internal states are 
necessary to describe the fluctuations of position of the motor, which are quantified 
by the so called randomness parameter introduced and measured in Ref. |48j . For 
this reason, a more refine model of kinesin should contain more than two internal 
states to account for the way the motor walks on the filament, which is by a succes- 
sion of binding and unbinding events of the two heads in a hand-over-hand fashion. 
To include that aspect, a model with 9 states - which can be reduced to 7 states 
for most cases - was proposed in Refs. [35j [49l [50] , where the 7 states describe 
the most significant chemical states of the two headed kinesin. The possible tran- 
sitions between these states can be represented by a network, which describes the 
mechano-chemical coupling in this motor. In this network representation, several 
cycles can be identified just like in the discrete model presented above. The model 
successfully accounts for many experimental results known for kinesins [51] . Using 
this framework, a diagram summarizing the thermodynamic modes of operation 
of the motor has been constructed in Ref. [55] . This diagram has similarities with 
our figure [5J but some differences are present due to the different role played by 
the mechanical and chemical cycles in the different models. 
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Many other works have used discrete or continuous stochastic models to 
analyze molecular motors : for instance a discrete model with 7 states has been 
developed for myosin V [S3]. A discrete model with only 3 states has been used 
to describe the fluctuations of position of nucleosomes along DNA in Ref. [SSJ . In 
Ref. [S3], the rotating motor iq-ATPase is described by a stochastic process for 
the angle of rotation of the motor, which is treated as a continuous variable and 
for the chemical states, which are treated as discrete states. More recently, the 
authors of this reference have developed a discrete version corresponding to their 
continuous model, which is a two states model very similar to the one discussed 
in this review [56) . 



3. Fluctuation relations in models of molecular motors 

Fluctuation relations quantify the exchanges of energy between a system and its 
environment when the system is in a non-equilibrium state |17H21j . These relations 
hold arbitrarily far from equilibrium in a regime where the usual thermodynamic 
laws - which hold only near equilibrium - do not apply. Since their discovery about 
a decade ago [TSJ [TBI HZ1 EH H] j there has been a growing interest to understand 
their importance and implications. One reason for the popularity of this topic has 
to do with the fact that these relations provide a fresh look at old fundamental 
questions, such as the origin of irreversibility or the second law of thermodynamics. 

For small systems (for which the fluctuations are large, in the sense that 
their magnitude can be of the same order as the average value), the fluctuation 
relations impose new constraints which go beyond the usual description of sta- 
tistical fluctuations. Many fluctuation relations have been verified experimentally 
using biopolymers, in particular the Jarzynski's relation [57], the Crooks relation 
[58) and the Hatano-Sasa relation [5 9) . Complementary experimental verifications 
have been carried out with colloidal particles in optical traps [50] HI] ■ Recently, 
a modified Fluctuation-Dissipation relation, related to the Hatano-Sasa relation 
has been verified for a colloidal particle in a nonequilibrium steady state [62] . 
All these experiments represent remarkable achievements, which confirm the va- 
lidity of the general framework of fluctuation relations in various experimental 
conditions. However, it may be worth pointing out that in all these examples, the 
experiments have been designed in order to verify the fluctuation relations. On 
the contrary, the case of molecular motors is particularly interesting since this is 
a system which was not designed for that objective. Molecular motors operate in 
a regime far from equilibrium, with fluctuation relations in some sense built-in in 
their natural mode of operation. For this reason, it is more appropriate to think 
about the fluctuations relations as thermodynamic constraints on the operation of 
the motors, which presumably theoretical models of molecular motors should obey 
[331 [32] . That of course would assume that the fluctuations relations are obeyed 
exactly in experiments. To our knowledge, quantitative experimental tests of fluc- 
tuation relations have not been carried out on molecular motors yet, although 
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there is some indication that, for instance, the data of |51j on single molecule 
experiments with kinesin is in agreement with the fluctuations relations. 

The dynamics of a molecular motor breaks the detailed balance condition, 
and leads to a non-equilibrium steady state characterized by the presence of non- 
zero currents, which are independent of time. For each current, one can associate a 
cycle, also called an irreversible loop. The construction of these cycles and the way 
they can be associated with currents is explained by a general theory for systems 
governed by a master equation [63 , 65 , 64 . One central result of this theory is the 
following relation 

n+(£) _ j+ _ A/kBT . . 

n-(£)~j-~ e ' (iZ) 

where II + (£) denotes the product of reaction rates associated with the different 
transitions within the cycle C in the clockwise direction, whereas II - (£) denotes 
the same product in the counter clockwise direction. We have denoted J + the 
number of cycles undergone by the motor per unit time in the clockwise direction, 
and J~ the number of cycles undergone per unit time in the counter clockwise 
direction, so that overall the cycle flux is J = J + — J~ . The quantity A is called 
affinity or thermodynamic force. This affinity can be defined as the derivative of 
an effective potential, experienced by a biased random walker that would exhibit 
the same dynamics 37 . When the detailed balance condition is satisfied, II + (£) = 
n _ (£), the effective potential is flat and A = J = 0. 




Figure 9. Cycles associated with the evolution of the motor in 
the discrete two-states model. Left: the cycle for the position vari- 
able n; the length n run by the motor corresponds to half the 
number of turns run around the cycle (the factor 1/2 has to do 
with the period of the motor which is twice the unit length of 
the lattice on which the motor evolves). Right: the cycle for the 
chemical variable y. 



In the simple case where the model contains a single cycle with only two 
states as in figure El the relation (fT2j) leads to the affinity A/ksT = —2^ where 
^> is defined by 
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and the corresponding current J is the average motor velocity (see Eq. I10[) : 

v — 2 6 — frufr b 

0? a + 73 b + tj a + tJ b ' 

In order to describe more precisely the dynamics of this system, let us consider 
Pj(n,t), the probability that the motor at time t is on the site i (= a, b) and at 
the position n (with x = nd where d is distance between sites a and b). This 
probability can be obtained for instance by integrating over the variable y in the 
quantity P n (y,t), which satisfies the more general master equation of Eq. (J9j) . 
Because of the periodicity of this problem, it is convenient to introduce generating 
functions Fi(\,t) = J2 n e _A "-Pi(n, t), which evolve according to : dtFi = MijFj, 
where .M[A] is the following 2x2 matrix constructed from the master equation 
satisfied by Pi(n,t): 



M[X] = 



-0j a —cj a e x tu~b + e X Lut, 



In the long time limit, the steady state properties of the motor can be obtained 
from the largest eigenvalue i?[A] of this matrix. Indeed when t — >• oo, 

(e- Xn )=J2m,t)~ex V ($t). (15) 

i 

The first derivative of $ with respect to A gives the average velocity v of the motor 
and the second derivative gives the diffusion coefficient of the motor. 

From the explicit expression of this eigenvalue, the following property may 
be derived 

0(A) = - A), (16) 
which is the Gallavotti-Cohen fluctuation theorem. Other equivalent forms of this 
relation can be obtained. One of them involves the large deviation function of the 
current v denoted G(v), defined in the long time limit, as: 

P{- = v) ~ e- G W* . (17) 

The analytical expression of this function, obtained in Ref. |37) , has a complicated 
non-linear expression in terms of the rates, but it satisfies a surprisingly simple 
relation: 

G(v)-G(-v) = *v. (18) 
This relation implies that the ratio of the probabilities to observe a velocity v or 
—v after a time t satisfies the relation: 

_ v l = e -* vt (191 

Using Eq. (|18l) . and the fact that G(v) and $(A) are Legendre transforms of each 
other, one recovers indeed the relation (|T6|) . 

The relations (fTBT) . (fT5]) and (|19l) are equivalent forms of a constraint imposed 
on the system by the Gallavotti-Cohen fluctuation theorem. This theorem itself 
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is a consequence of the time-reversal symmetry of the physical laws involved in 
this model. This symmetry is a fundamental property that does not depend into 
on the details of the system and therefore, in this sense, the fluctuation theorem 
appears as a universal requirement, just as thermodynamic constraints are univer- 
sal for systems at equilibrium regardless of their microscopic structure. Of course, 
universality does not imply that the constraints are easy to find and to formu- 
late explicitly for a given problem. Here, the relation (|19|) is an explicit prediction 
derived from the fluctuation theorem for molecular motors. It would be of great 
interest to verify this relation experimentally. Conversely, this equation can be 
used to measure the affinity ^P, and therefore to access indirectly the microscopic 
rates, for a given molecular motor. 

3.1. Modeling processivity at the single motor level 

Experiments on single molecular motors depend on an important property of these 
motors called processivity. Molecular motors like kinesins, which can hydrolyze a 
large number of ATP molecules before detaching from microtubules, are called 
processive, whereas those like myosins II, which detach and reattach frequently 
from actin filaments and are called non-processive. There are several ways to de- 
fine processivity, either it can be defined as the average lifetime of the motor on 
the filament, or from the average length spanned by the motor or from the average 
number of ATP molecules consumed before detaching. In single molecule experi- 
ments (66] [51] , the dependence of the run- length of a single kinesin on load and 
on the ATP concentration has been studied. On the theoretical side, the average 
lifetime of a molecular motor as function of load has been studied in .67[ using 
the flashing ratchet model. Here we focus on the run-length, for which we derive 
a simple expression within the discrete two states model presented above. Using a 
similar theoretical approach, the average time before observing a backward step in 
a discrete model of a molecular motor has been studied in Ref. [33] . More generally, 
in a network of discrete chemical states, there are well known methods to calculate 
the average lifetime of a random walker in the presence of absorbing states [65] . 
These methods can be used not only to calculate the lifetime of the motor on the 
filament, but also the dwell times associated with the motor steps as shown in 
Ref. H5]. 

The detachment of the motor from the filament can be represented by a 
detachment rate re, which depends on the local site on the filament visited by the 
motor. Since in the discrete model presented above, the state b is the high energy 
state and a a low energy state, we assume for simplicity that detachment only 
occurs from site b, which corresponds to the maxima of U\(x) in figure [3] Because 
of this detachment, the motor can be no longer only in states a or &, thus we need 
to modify the master equation in order to conserve probability at all times. This 
can be done by adding an extra state corresponding to the unbound motor state, 
which is an absorbing state. We define the generating functions Fi(X,t) as before, 
but now the matrix of evolution of these generating function is the following 3x3 
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matrix: 



M[X] 



K 



This matrix has three eigenvalues fi 2 and 0, which is associated with 
the absorbing state. The corresponding eigenvectors are IM2) and |c). If the 
initial state vector is |-F(A, 0)) = A\fj,i) + B\^i 2 ) + C\c), the state vector at time t 
is \F{\,t)) — Ae^ lt \^i) + Be^ 2t \^ 2 ) + C\c). Since /xi and \i 2 are strictly negative, 
at time t ->■ 00, (e~ An ) = (0|F(A,oo)) = C(A), with (0| = (1,1,1). Note that 
C(A) contains all the moments of the run length, and in particular the average run 
length, which is, in units of the lattice period, 

<n> = -^U=o. (20) 

The function C(A) can be calculated by projecting the left eigenvector associated 
the eigenvalue 0, (c|, onto the initial state vector \F(X, 0)). If the initial state vector 
is along the unit vector e x , then one obtains 

r(x) _ , + k« e ~ A m , 

{ ' ™ a + W^bjl - e 2A ) + U a tJ b (l - e-^) ■ 

If the initial state vector is along the unit vector e x , then 

(n = ■ h2 . (22 

UJ a + LOb KUJ a 

If the initial state vector is along the unit vector e y , then 

(n)=2 ^ 6 ~ t7at7b . (23) 

KU) a 

Thus the part of the average length which is independent on the initial condition 
is 

<») = 2 ^" t7at7t = 4-, (24) 

where v is the average motor velocity defined before and Pb the stationary proba- 
bility to be in state b when k = 0. 

Note that it has been assumed implicitly that the motor runs in the positive 
direction so that by construction v > 0, (n) is positive and has the familiar form 
obtained above. The application of fluctuation relations to non-processive motors 
has not been discussed in the literature to our knowledge. We believe that a fluc- 
tuation relation will be obeyed only if a reattachment process is taken into account 
in the model. In this case, there is no absorbing state anymore. 
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3.2. Mechanochemical coupling for the discrete model 

We have so far only discussed the form of fluctuation relations for models contain- 
ing a single cycle. It is well known that models with at least two cycles must be 
introduced to account for experimental data such as those represented in Figure [SJ 
In order to discuss more general fluctuation relations, and to compute in a simple 
way the chemical current, r, associated with the average ATP consumption rate, 
it is necessary to include in the description of the state of the motor, a chemi- 
cal variable y associated with the average number of ATP consumed as done in 
Eq. (j9)). With the notations, u> l a — ~uJ l a + tu L a and ui a — w a + to a (and similarly for 
site 6), one obtains the chemical current 

r = -2-5 b - -, (25) 

in agreement with the cycle representation of figure |H] and with the formula (|11[) 
above. 

When the form of these rates is explicitly given in terms of the normalized 
force / and the normalized chemical potential A/x, a new reformulation of Eq. (|12[) 
is obtained: 

fc B Tln = F e (2d) -Afi(l- , (26) 

with I, I' = 0, 1. This equation can be understood as a statement of the first law 
of thermodynamics at the level of elementary transitions [35] [28]. Indeed, it is 
possible to associate the left hand side of this equation with the heat released by 
the motor into the environment (treated as a reservoir at the same temperature) 
during transitions (Z, /'). The right hand side can then be interpreted as the differ- 
ence between the mechanical work —F e (2d) and the variation of chemical energy 
Ap, (I — I') for these transitions. The variation of internal energy is zero in this case 
since only transitions involved in a cycle are considered. Note that similarly to this, 
the generalized detailed balance condition of Eq. ([8]) can be also interpreted as a 
statement of the first law at the level of elementary transitions. 

Following the same steps that lead to the fluctuation relations for models 
with one cycle but now for a model with two cycles, one arrives at the following 
relation, similar to Eq. (1 1 G[) : 



0(A,7) = 0(-*-A,-x-7). ( 27 ) 
with new affinities — \£ and — x associated with the mechanical and chemical cycles. 
These affinities represent a part of the expression of the entropy production rate 
of the motor [37], which also satisfies a fluctuation relation different from that of 
the currents but very much related to it. 

We emphasize that the relation (|2T|) involves both the mechanical and the 
chemical activities and that, here, a symmetry relation of the type of Eq. ([16]) for 
the mechanical variable alone is not satisfied. It is therefore essential, in deriving 
fluctuation relations, to take into account and include all internal variables that 
are coupled with one another under time-reversal. Leaving aside some relevant 
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degrees of freedom would manifest itself as an apparent violation of this symmetry 
and would wrongly be interpreted as a breakdown of the fluctuation theorem. 

3.3. Flashing ratchet model on a continuous space 
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Figure 10. Normalized eigenvalue T<d(fr)) of the flashing ratchet 
model as a function of r\ = A// (with r = a 2 /Do), for a nor- 
malized force / = 5 (top two curves) and / — 10 (bottom two 
curves). The solid curves correspond to the case where the tran- 
sition rates between the internal states satisfy detailed balance, 
which leads to the Gallavotti-Cohen symmetry, i.e. to the sym- 
metry with respect to rj — —1/2. The curves with filled symbols 
(/ = 5) and empty symbols (/ = 10) correspond to the case 
where the detailed balance is broken with constant transition rates 
uii(x) = ui2{x) = 10r _1 and the same potentials. 

We show in this section how the ideas developed in the previous section for the 
discrete two states model can now be extended to more general continuous models, 
such as the flashing ratchet. There are several reasons for considering continuous 
models as a substitute for discrete models: first of all, continuous models contain all 
the possible discrete models as limiting cases, second, the way to describe the effect 
of force on the motor is unambiguous for continuous models, and third, there are 
effects such as fluctuations which are not always well captured by discrete models. 

We provide in this section an analytical proof that the flashing ratchet obeys 
a Gallavotti-Cohen symmetry [68], using a technique inspired by 17, 22J. We also 
analyze numerically this point by calculating the eigenvalue associated with the 
evolution matrix of the generating functions of the currents. 

3.3.1. The purely mechanical ratchet. Before considering the case of the flashing 
ratchet with a mechanical and a chemical variable, it is helpful to look first at 
a purely mechanical ratchet, which has been used to describe in particular the 
translocation of a polymer through a pore [69] . In this model, one considers a 
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random walker in a periodic potential subject to an external force F (model I) 
[701 I3"9"] . The corresponding Fokker-Planck equation is 

~dP U'(x)-F 



dP_ _ d_ 

dt dx 



dx knT 



-P 



(28) 



where U (x) is a periodic potential U(x + a) = U (x) and a is the period. This 
equation describes the stochastic dynamics of a particle in the effective potential 
U e ff(x) — U(x) — Fx. By solving Eq. (|28p with periodic boundary conditions 
[501 170] , it can be readily proven that the system reaches a stationary state with 
a uniform current J in the long time limit. This current is non- vanishing if a non 
zero force is applied. When F — 0, there is no tilt in the potential, J = and the 
stationary probability is given by the equilibrium Boltzmann-Gibbs factor. 
Similarly to the discrete case, we introduce the generating function 

Fx((, *) = E exp (A(C + n)) P((n + C)o, i) . (29) 

a 

The time evolution of this generating function F\ is obtained by summing over 
Eq. (j2"5)l . This leads to the following equation: 

d -E^A = C {\)F x {(:,t) 1 (so) 

where the deformed differential operator C{\) acts on a periodic function i) 
of period 1 as follows: 

a 2 <9 2 <I> 8 /- \ 5$ 

= |£ + ^ (t/*„/*) - 2A^ - AE^,* + A 2 $ , (31) 

where U' e ft = ad x ll e f / /ksT and the left hand side of Eq. (I3T1) is proportional to 
the inverse of the characteristic time r = a 2 /Do. 

The operator £(A) has the following fundamental conjugation property: 

:bT £{\) f e -u(x)/k B T^ = £t ( _j _ A) $ (32) 

This property implies that the operators C(X) and (— / — A) are adjoint to each 
other, and thus have the same spectrum. If we call 0(A) the largest eigenvalue of 
£(A), we obtain from Eq. ([32| that 6(A) satisfies the Gallavotti-Cohen symmetry: 

9(A) = 9(-/ - A). (33) 

In fact, this symmetry holds for all eigenvalues. For the special case / = 0, the 
conjugation relation ([32| reduces to the detailed balance property [22]. One can 
note that this proof of the Gallavotti-Cohen symmetry does not require explicit 
knowledge of 0(A). In the discrete minimal ratchet model, an explicit analytical 
expression could be obtained for this quantity. In the continuous case, this is no 
longer the case but O(A) can be calculated numerically. We have done this by first 
discretizing the operator £(A) and then calculating its largest eigenvalue using the 
Ritz variational method 681. A similar method has been used in Ref. [7T] for the 
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cosine potential. We note that our numerical approach can handle any form of 
potential. 

3.3.2. The flashing ratchet. We now present the extension of the Gallavotti-Cohen 
symmetry to the case of the flashing ratchet model, which should include both 
the mechanical and chemical currents [3H1 [2E] ■ When the switching rates satisfy 
a detailed balance condition, which is for instance the case when A/x = 0, the 
symmetry is indeed present as shown in the solid curves of Fig. 1101 In the general 
case however, where the normalized force / and chemical potential A/i are both 
non-zero, the relation (J7J is no longer satisfied and the Gallavotti-Cohen relation 
(pi3")) is not valid. This is shown in the curves with symbols in Fig. [TU] where for 
simplicity we took constant switching rates u)\ — u>2 — 10r _1 . For all the curves 
of this figure, a sawtooth potential U%, and a potential Ui constant in space have 
been chosen. The breaking of the symmetry of Eq. (|33[) can be interpreted as a 
result of the existence of internal degrees of freedom. Although other mechanisms 
exist which lead to violations of fluctuations relations as discussed in Ref. [71] , this 
case appears to be rather generic. 
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Figure 11. For the model described by Eqs. (|34ll35j) . the nor- 
malized eigenvalue tA(/?7, — A/i/2) is shown as function of rj. The 
dashed curve corresponds to / = 5 and A/i = 0, the solid curve 
corresponds to / = 5 and A/i = 10, and the dotted curve corre- 
sponds to / = 2 and A/i = 10. The symmetry is recovered in all 
cases in this description which includes both the mechanical and 
chemical degrees of freedom. 

Let us now introduce the probability density Pi(x,q;t) associated with the 
probability that at time t the ratchet is in the internal state i, at position x and 
that q chemical units of ATP have been consumed. The evolution equations for 
this probability density is obtained by modifying Eqs. ([3]) after taking into account 
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the dynamics of the discrete variable q. We have 
dPi(x,q,t) 



dt 

dP 2 (x,g,t) 
dt 



= (A -u)i(x))Pi(x,q,t) (34) 

+ lo^{x)P 2 [x,q + l,t) + uj 2 ) (x)P 2 (x, q, t) 

= {L 2 -U2{x))P 2 {x,q,t) (35) 

+ ul{x)P x [x, q, t) + u\ (x)P 1 {x,q-l,t). 

We use a notation similar to that of Ref. 37 , where oj\(x) denotes the transition 
rate at position x from the internal state i with l(= —1,0,1) ATP molecules 
consumed. This leads to 

(36) 
(37) 
(38) 
(39) 

We also have lo\{x) = tOi{x) + u{(x) and ^2(2) = w§(a;) + aj 2 _1 (x). The operators 
C\ and £2 act on a function 3> as 

As above, we introduce two generating functions i*i,A,7 and F2,\,j, depending 
on two parameters A and 7 which are conjugate variables to the position of the 
ratchet and to the ATP counter q. We have for i = 1, 2, 

JWC *) = E e7? E e A(C+n) ^(a(C + n), 9; i). (41) 

g n 

The evolution equation for these generating functions is obtained from Eq. (|35[) as 



W °(z) 


= 


-fx)/k B T 


u° 2 (x) 




-fx)/k B T 


oj{(x) 


= ^ 


-fx)/k B T+Afi 


' 2 - 1 W 




-fx)/k B T _ 



a <vft,x.w- £(A ' 7> u,;::j- (42) 

with the operator £(A, 7) decomposed as 

£(A, 7 )=P(A)+AA( 7 ), (43) 

with T>(\) the diagonal matrix diag(£i(A) — ui\, £2 (A) — w 2 ), where the deformed 
operators C\ (A) and £2 (A) have the form written in Eq. (|31[) with ?7 e j j (x) given 
by Ui(x) — Fx for i = 1,2, respectively. The operator A/"(7) is defined as 

"° + f e " T )- (44) 

Consider now the diagonal matrix Q defined by diag(e _l/l /' CBT , e~ U2 ^ k]BT ). By di- 
rect calculation, one can verify that Q~ 1 N{^)Q = A/"' (—A// — 7) . From Eq. 
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one obtains Q l T>{i)Q = (— A/i — 7). By combining these two equations, we 
conclude that 



where A(A, 7) is the largest eigenvalue of £(A, 7). This relation is the equivalent 
of Eq. P7|) . which was derived for the discrete model. If we consider only the 
mechanical displacement of the ratchet, the relevant eigenvalue 9(A) is given by 
0(A) = A(A,0), which clearly does not satisfy the fluctuation relation of the form 
Eqs. (|16H19p as shown in Fig. [TO! In Fig. [TTJ we have computed A(fr), — A/i/2) for 
the same potentials and with rates oj\(x) of the form given above with oj(x) = 5r _1 
and (j>{x) = 10t _1 . We have verified that in all cases the symmetry of Eq. pp]) 
holds. 

To conclude this section on the continuous flashing ratchet, we emphasize 
the following two points: (i) we have proved that the flashing ratchet satisfies 
the fluctuation theorem without having to adjust any parameter in the system to 
enforce the validity of this theorem. The only constraints on the switching rates 
were given a priori from thermodynamics and kinetic theory and these require- 
ments are always taken into account in the very definition of the model (see e.g. 
[4Ti]). What we have shown is that these thermo- kinetic constraints are enough to 
imply the Gallavotti-Cohen symmetry, which itself has far reaching consequences 
on the model, (ii) In order to derive the fluctuation theorem, all relevant micro- 
scopic degrees of freedom must be involved. For example, in the flashing rachet 
model, the position variable alone does not obey the fluctuation theorem and the 
chemical variable that counts how many ATP molecules have been consumed by 
the motor during its displacement has to be taken into account. The Gallavotti- 
Cohen symmetry thus leads to global mechano-chemical constraints on the modes 
of operation of the motor. 



4. Conclusions 

A first simple and useful message to take from this study is that the dynamics 
of a molecular motor can be described by the evolution of a random walker in 
an effective potential U e f f{x,y) where x is the mechanical variable and y is the 
chemical variable [75j . The periodicity of the potential along x and y implies that 
the potential has an egg-carton shape. 

The symmetry of the fluctuation relations for the currents is valid in general 
for the flashing ratchet model only when internal degrees of freedom arc taken into 
account. This raises a fundamental question concerning the validity of fluctuations 
relations and their applicability to other types of ratchet models [751 [73J HI |3H] . 
More generally, other mechanisms exist which are known to produce deviations 




(45) 



(46) 



Fluctuation Relations for Molecular Motors 



25 



from fluctuations relations [74 , and it would be valuable to know whether fluc- 
tuations relations can always be restored by enlarging the phase space and by 
modifying the dynamics accordingly. 

On the experimental side, it would be very interesting to investigate fluctu- 
ations relations for molecular motors using single molecule experiments, in a way 
similar to what was achieved in colloidal beads or biopolymers experiments [57j - 
|62j . Using fluorescently labeled ATP molecules, recent experiments with myosin 
5a and with the i<i-ATPase rotary motor, aim at simultaneous recording of the 
turnover of single fluorescent ATP molecules and the resulting mechanical steps 
of the molecular motor [76] . These exciting results indicate that a simultaneous 
measurement of the values of the mechanical and chemical variables of the motor 
is achievable in practice, and therefore from the statistics of such measurements 
it may be possible to obtain the distribution of probability to find the motor at a 
specific position and with a specific number of molecules of ATP consumed. With 
enough statistics, one could thus in principle verify Eq. (|27p. Such an experimental 
verification would confirm that the Gallavotti-Cohen symmetry is a fundamental 
constraint that plays an essential role in the mechano-chemical coupling of molec- 
ular motors. 

Finally, besides the Gallavotti-Cohen fluctuation theorem, many exact non- 
equilibrium relations have been discovered during the last decade, the most famous 
one being the Jarzynski identity |16j and its generalization by Crooks |19) . These 
identities, originally derived for systems being driven out of a state of thermo- 
dynamic equilibrium, have been extended by Hatano-Sasa to systems prepared 
in non-equilibrium stationary states and following Markovian dynamics [77 . This 
more general fluctuation relation leads with a proper choice of observables to gen- 
eralizations of the well-known fluctuation-dissipation known for systems close to 
equilibrium [78, 79, 80,;81]. These generalized fluctuation-response relations hold 
for systems prepared in non-equilibrium stationary states and following Markovian 
dynamics. The implications of the Crooks fluctuation theorem for kinesin has been 
analyzed in Ref. |82) , while the implications of generalized fluctuation-response re- 
lations for molecular motors have been considered recently J53J HI] • 

All these various relations can be interpreted as universal constraints that 
have to be obeyed by systems far from equilibrium regardless of their detailed 
structure. It would be of great interest to explore the consequences of these rela- 
tions in the field of ratchet models (both at the single motor level and at the level 
of many motors) and to draw from them some measurable predictions that could 
be verified experimentally on molecular motors. 
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